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ABSTRACT 

We study the effect of dissipational gas physics on the vertical heating and thick- 
ening of disc galaxies during minor mergers. We produce a suite of minor merger 
simulations for Milky Way-like galaxies. This suite consists of coUisionless simulations 
as well as hydrodynamical runs including a gaseous component in the galactic disc. 
We find that in dissipationless simulations minor mergers cause the scale height of 
the disc to increase by up to a factor of approximately two. When the presence of 
gas in the disc is taken into account this thickening is reduced by 25% (50%) for an 
initial disc gas fraction of 20% (40%), leading to a final scale height Zq between 0.6 
and 0.7 kpc (for a sech^ profile), regardless of the initial scale height. We argue that 
the presence of gas reduces disc heating via two mechanisms: absorption of kinetic 
impact energy by the gas and/or formation of a new thin stellar disc that can cause 
heated stars to recontract towards the disc plane. We show that in our simulations 
most of the gas is consumed during the merger and thus the regrowth of a new thin 
disc has a negligible impact on the zq of the post merger galaxy. Final disc scale heights 
found in our simulations are in good agreement with studies of the vertical structure 
of spiral galaxies where the majority of the systems are found to have scale heights of 
0.4 kpc ^ Zq < 0.8 kpc. We also found no tension between recent measurements of the 
scale height of the Milky Way thin disc and results coming from our hydrodynamical 
simulations. Even if the Milky Way did experience a recent 1:10 merger it is possible 
to reproduce the observed thin disc scale height, assuming that the disc contained 
at least 20% gas (similar to the gas fraction today) at the time of the merger. We 
conclude that the existence of a thin disc in the Milky Way and in external galaxies 
is not in obvious conflict with the predictions of the CDM model. 

Key words: Galaxy: disc, evolution, structure - galaxies: disc, evolution, interactions, 
structure - methods: numerical, N-body simulation 



1 INTRODUCTION 

The cold dark matter (CDM) theory provides a success- 
ful framework for understanding structure formation in 
our Universe. Within this paradigm, dark matter first col- 
lapses in small haloes, wh ich merge to form progressively 
larger haloes over time (e.g. lWhite fc Rees|[l97 8l: lDavis et al.l 
[l98a). Major (near-equal mass) and minor (unequal mass) 
mergers are a generic feature of structure assembly in the 
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CDM or hierarchical picture, and numerical simulations 
have shown that many of these merging structures sur- 
vive within the vir i alized reg ion of larger host haloes (o^. 
Klvpin et aP Il999l: iReed et al . 2005; Diomand ct al. 2008 



I- ■ - . . .^^^^^uji — — ; 1 ~ --^j 

20081 ). These merger events are now widely 



believed to be responsible for shaping many galaxy proper- 
ties. Major mergers play an important role in tr ansforming 
disc-dominated spira l galaxies into spheroids (e.g. lHernquistI 
1 19931 : iNaab fc Burkert .2003,) and triggering episodes of 
enhance d star formation (SF) and a c tive galactic nuclei 
(AGN) jHernauistl [19891 : iBarnej Il992l : iMihos fc Hernguistl 
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1994 iBarnes fc Hernguis^ Il996l: iMihos fc Hernguisd Il996l : 



On the observational side, there is a general consen- 



Springel et all l2005l : [Johansson et alj |2009| ). Minor merg- sus that ~70% of Milky- Way-mass 10^^ Mr) haloes host 

a disc dominated, la t e type galaxy |Weinmann et al. 200' 



ers may explain the origi n of th ick discs (Abadi et al 
20031 : iBrook et al.l i2004: R ead et alJ lioOS: Kazantzidis ct al. 
20081 ). and the diffuse stellar halo around galaxies may be 
produced via tidal destruction of merging sat e llites (e.g 



BuUock fc JohnstonI l2005l : iMurante etal] I2OO4I : iBeU et all 
20071 ). 



On the other hand, this large population of merging 
satellites has raised the question of whether mergers are 
too common in the CDM scenario. Some studies have ques- 
tioned whether CDM models can produce a large enough 
population of "bulgeless" discs (or syst ems with very low 
bulge-to-total ratios B/T < 0.2) (e.g. iGraham fc WorlevI 
l2008l : IWeinzirl et al.l |2009| ) and whether thin, dynamically 
fragile discs such as the one observed in the Milky Way 
can survive thi s bom bardment by incoming satellites (e.g . 
Toth fc Ostrik cr 1992: Quinn et al.lll993l : IWalker et al]|l99d : 



WyselbOOlflBournaud et al.ll2007n . Clearly, there are three 



main aspects to settling this question: first, we must under- 
stand the statistics of galaxy mass accretion and merger his- 
tories in a CDM universe; second, we need to understand the 
physics of how galaxies are transformed by mergers with var- 
ious mass ratios, orbits, and other parameters; and finally, 
we need accurate and unbiased statis tics for the observed 
popu lations of the relevant objects (cf. iKormendv fc Fisher! 
l2008h . 

The availability of large, high-resolution N-body sim- 
ulations has made pro gress possible on the first of these 
aspects. For example, Istewart et al.l l|2008l ) recently car- 
ried out a detailed study of statistics of merger events in 
galaxy-sized dark matter haloes. They confirmed previous 
results based on semi-analytic methods l|Purcell et al.|[2007l : 
IZentnei:j|2007l ). finding that the majority of the mass delivery 
into a dark matter halo of mass Mh is due to systems with 
masses Mgat = (0.03 — 0.3)Mh. They also found that a large 
fraction (95%) of Milky- Way sized haloes have accreted a 
satellite with a virial mass comparable with the total mass 
of the Milky Way disc and approximately 70% have accreted 
an object with more than twice this mass since z ~ 1. 

There have been numerous studies based on collision- 
less simulations that have tried to quantify the effects of 
these minor m ergers on the thickness and stability of stel- 



lar discs (e.g . Velazquez fc White 1999; Fon t et all I2OOII: 



Benson et al.ll2004: Gauthier et al., ,2006.: Kazantzidis et al] 



20081 . l2009l : IVillalobos fc Helmill2008l) . This work has demon- 



strated that the answer depends quite sensitively on the 
mass ratio of the satellite to the primary; there seems to be a 
consensus that the main danger to thin discs is not from the 
ubiquitous mergers with very small mass ratios (less than 
1:10) but rather from the rarer, yet still fre quent, events with 
larger mass ratio (^1:10). For example [Kazantzidis et al] 
l|2008h studied the effects of mergers with (0.2 — l)Mdisc 
and concluded that thin discs could s urvive such m e rgers . 
Employing dissipationless simulations, iPurcell et al.l (|2009l . 
P09) studied the response of a fully formed Milky Way-like 
stellar disc within a ~ 10^^ M0 DM halo to mergers involv- 
ing satellites with a total mass Maat ~ 10^^ Mq ~ 3Mdisc- 
They came to the conclusion that regardless of the orbital 
configuration of the merger these events transform the discs 
into structures that are roughly three times thicker than the 
observed thin disc component of the Milky Way. 



a disc dommated, la t e type ga 
van den Bosch et all l2007l : jc 



Choi et al. 2007; Park ct al 



20071 ). While it is also well known that the majority of the 



mass in the disc of our Milky Way resides in a thin com- 
ponent (e.g. lJuric et al.l l2008:l. how typical this situation is 
for other disc galaxies in the Milky Way's mass range is less 
certain. Quantifying the vertical thickness of the discs of 
external galaxies has been attempted in a few studies (e 
[Schwarzkopf fc Dettnia3l2000l : lYoachim fc DalcantonllioQl 
but is challenging because of small sample sizes, inclination 
effects and extinction. 

P09 suggested two possible explanations for the exis- 
tence of a thin disc in our Galaxy: one possibility is that the 
Milky Way is not a representative case and that perhaps it 
has had an unusually quiet accretion history for a halo of 
its mass. A second explanation is related to the fact that all 
the numerical studies performed so far have only considered 
the dissipationless components in the galaxy (dark matter 
and stars), neglecting the presence of a dissipative gas com- 
ponent in the disc. However, the inclusion of gas physics is 
known to play an important role in stabilizing galactic discs. 
Numerical simulations have shown that gas is important 
for the survivabilit y and t he re ^rowth of stellar discs dur- 
ing major merg ers (Barnes 2002; Springcl & Hcrnquist 200f 
Robertson et al. 2006; Naa b ct al.. .2006 ; Scannapicco et al 



20091 : iGovernato et al.ll2009l ), and that the presence of gas 
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disc-type merger progenitors greatly suppres ses the forma- 
tion of a post-merger spheroidal component iHopkins et al.l 
(2009allB). It is reasonable to expect that the presence of gas 
in the progenitor disc could also have an impact on the ef- 
ficiency of the disc heating and thickening: gas may be able 
to absorb some of the kinetic impact energy of the merging 
satellite, and then radiate this energy away by cooling; or, 
gas may be able to cool and reform a new thin disc after the 
merger, forcing heated stars to contract again onto the disc 
plane. 

The main goal of this paper is to study in detail how the 
presence of a dissipational gas component affects disc thick- 
ening in minor mergers. For the first time, we address this 
problem using a suite of high-resolution, fully hydrodynam- 
ical numerical merger simulations, run with the Smoothed 
Particle Hydrodynamics (SPH) code GADGET-2. We con- 
sider simulations with and without gas, with a variety of 
initial orbital parameters, and with different values of the 
star formation efficiency parameters. We then re-examine 
the issue of whether the cosmologically expected rate of mi- 
nor mergers presents a problem for CDM in light of our new 
results and the available observations. 

The remainder of the paper is organized as follows: in 
Section [2. II we provide a brief summary of the GADGET-2 
code and describe our simulations as well as the initial con- 
ditions for our progenitor galaxies. In Section [3] we present 
our main results, focusing on the difference between the dis- 
sipational and dissipationless simulations as well as on the 
inffuence of different initial conditions and star formation 
parameters on the final thickness of the disc. Finally in Sec- 
tion|3]we summarize and discuss our results, compare them 
with observations, and present our conclusions regarding the 
issue of thin disc survival in a CDM universe. 
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Table 1. Parameters kept constant for all simulations. Masses are in units of 10^'' Mq; scale and softening 
lengths are in units of kpc and pc respectively. 



System 


Afvir 






rd 


fb 


c 




-^disc 


-^bulgc 


EDM 


^disc 


Primary 


100 


2.4 


0.600 


3.0 


0.5 


9.65 


4 000 000 


1 000 000 


500 000 


100 


50 


Sat 


10 


0.0 


0.063 


0.0 


0.3 


11.98 


900 000 





100 000 


100 


50 



2 NUMERICAL SIMULATIONS 



2.1 Numerical Code 

We make use of the parallel TreeSPH-code GADGET-2 
l|Springelll2005l ) in this work. The code uses Smoothed Parti- 
cle H ydrodynamics (SPH: ILucvII 19771 : [Cingold fc MonaghanI 
1 1977 : Monaghan 19 92,) to evolve the sas using an entropy 
conserving scheme (jSpringel fc Hernguistl [20021 ). Radiative 
cooling is implemented f or a primor d ial m ixture of hydro- 
gen and helium following iKatz et al.l l|l996l ) and a spatially 
uniform time-indepen dent local U V background in the op- 
tically thin limit (Ha ardt fc Mad au 1996) is included. 

The SPH properties of the gas particles are averaged 
over the standard GADGET-2 kernel using ~ 64 SPH par- 
ticles. Additionally the minimum SPH smoothing length is 
required to be equal to the gravitational softening length in 
order to preven t artificial stabilization of small gas clumps at 
low resolution (|Bate fc Burkertll 19971 ). All simulations have 
been performed with a high force accuracy of Oforcc ~ 0.005 
and a tim e integration a ccuracy of j^acc = 0.02 (for further 
details see lSpringeill2005l ). 

Star formation and the associated heating by super- 
novae (SN) is modelled follo wing the sub-resolution mult i- 
phase ISM model described in lSpringel fc Hernguistl (|2003l ). 
The ISM in the model is treated as a two-phase medium with 
cold clouds embedded in a hot component at pressure equi- 
librium. Cold clouds form stars in dense (p > pth) regions on 
a timescale chosen to match observations (jKennicutt 1998) . 
The threshold density pth is determined self-consistently by 
demanding that the equation of state (EOS) is continuous 
at the onset of star formation. We do not include SN-driven 
galactic winds nor feedback from accreting black holes (AGN 
feedback) in our simulations. 

In our "fiducial" runs, we adopt the standard parame- 
ters for the multiphase feedba ck model in order to match th e 
Kennicutt Law as specified in ISpringel fc Hernguistl (|2003h . 
The star formation timescale is set to =2.1 Gyr, the cloud 
evaporation parameter to Ao = 1000 and the SN "temper- 
ature" to TsN ~ 10* K. In order to test whether the val- 
ues of these (uncertain) parameters affect our results, we 
reran two simulations with parameters that are a factor of 
4 larger (t'j = 8.4 G yr, Ao = 4000, Tsn = 4 x 10* K). As 
ISpringel et all (|2005l ) noted, this choice of parameters re- 
sults in a star formation rate (SFR) of ~ 1 Mq yr~^ for a 
Milky Way-like galaxy and gives better agreement with the 
long gas consumption timescale inferred for the Milky Way. 
However, these parameters yield a SFR that lies slightly be- 
low the Kennicutt Law. 



2.2 Galaxy Models 

We apply the method described in ISpringel et al.l l|2005l ) to 
construct both the central galaxy (called the primary here- 
after) and the merging satellite systems. Each primary sys- 
tem consists of gas and stellar discs with radial profiles de- 
scribed by an exponential, with masses Mgas and Mdisc , 
and a spherical bulge of mass Mh embed ded in a dark mat - 
ter halo of mass A'Uir ■ The halo has a iHernguisO (|l990l ) 
profile with a scale radiu s a corre sponding to a Navarro- 
Frenk- White halo (NFW; iNavarro et a l. 1997) with a scale 
length of Ta and a concentration para meter c — r^ii/ra. We 
use the results of iMaccio et al.l (|2008| ) to compute halo con- 
centration as a function of virial mass. 

The scale lengths of the exponential gaseous and stel- 
lar discs are assumed t o be eq ual, and are determined using 
the model of IMo et all l|l998l ), assuming that the fractional 
angular momentum of the total disc ja ~ (Jgas + Jdisc)/J'vir 
is equal to the global disc mass fraction = (Afgas -I- 
A/disc)/A/vir for a constant halo spin A . This is equivalent 
to assuming that the specific angular momentum of the ma- 
terial that forms the disc is the same as that of the initial 
dark matter halo, and is conserved during the process of disc 
formation. 

The vertical structure of the stellar disc is described by 
a radially independent sech^ profile with a scale height zq, 
and the vertical velocity dispersion is set equal to the radial 
velocity dispersion. The vertical structure of the gaseous disc 
is computed self-consistently as a function of the surface den- 
sity by requiring a balance of the galactic potential and the 
pressure given by th e EOS . The stellar bulge is constructed 
using the iHernguistI l|l990l) profile with a scale length rb . 

Satellite systems are constructed using the same 
method adopted for primary systems but they consist only 
of a dark matter halo and a stellar bulge; neither stellar nor 
gaseous discs are included. 



2.3 Simulation Parameters 

We construct a set of primary systems, each with a virial 
mass of A/vir = IO^^A/q containing a disc and a bulge, and 
a satellite system with a virial mass of A/vir = 10" A/0 con- 
taining only a bulge. Since we wish to study a typical galaxy 
in a A/vir = IO^^A/q halo, rather than use the specific pa- 
rameters of the Milky Way, we use th e average stellar-to-halo 
mass ratio derived by iMoster et al.l (|2009l ) to compute the 
stellar m ass of eve r y syst em (both primaries and satellites). 
The Moster et all l|2009l ) constraints were derived empiri- 
cally by asking how the population of dark matter haloes 
and sub-haloes predicted by CDM must be populated with 
galaxies in order to reproduce the observed galaxy stellar 
mass function, and are in excellent agreement with con- 
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Table 2. Parameters for the different simulation runs. Masses are in units of IO^^Mq, and scale height and 
softening lengths are in units of kpc and pc respectively. The star formation timescale (t^ ) is expressed in 
Gyr. The first three entries are for isolated galaxies, and the remainder are for mergers. 
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U 
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1 \y 1 n8 r/' 
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0.2 
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125 000 


70 


0.034 


0.40 




2.1 
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1 X 10*K 


IC 


0.4 


1.6 


4.0 
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70 


0.036 


0.40 
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2.1 
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1 X lO^K 


MA60 


0.0 
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2.4 








0.033 


0.40 


60° 
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1 X 10*K 


MB60 
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0.6 
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0.40 


60° 
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1 X 10*K 
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MA60T 
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60° 


8.4 


4000 


4 X 10*K 



straints from other methods such as galaxy clustering, satel- 
lite kinematics and weak lensing. Note that the stellar mass 
we get for a typical galaxy in a Mvir = W^^Mq halo is sig- 
nificantl y less than what is found for the Milky Way and 
M31 (cf. lKlypin et al.ll2002l ). 

All primary systems have a concentration parameter 
of c = 9.65 and a stellar mass of M.,pri = 3 X 10^° Mq. 
Distributing 80% of this stellar mass into the exponential 
disc yields a stellar disc mass of Mdiac = 2.4 X 1O^°M0 and 
a bulge mass of Mb,pri = 6 x 10® M©. We assume a bulge 
scale length of rb,pri = 0.5 kpc. 

For models that also include a gaseous disc compo- 
nent we add Mgas,20% = 0.6 x 10^"Mq and Mgas,4o% ~ 
1.6 X 10^" Af© such that t he gas fraction in the disc is 20% and 
40% in the two cases (cf . IStewart et al]|2009l ) . The total disc 
mass is then Md,o% = 2.4x 10^°Mq, Md,2o% = S.Ox lO^^M© 
and Md,4o% = 4-0 x IO^^Mq for a gas fraction of 0%, 20% 
and 40% respectively. We decided to keep the total stellar 
mass constant in order to have a more direct comparison be- 
tween the dissipational and dissipationless cases. We fix the 
disc scale ra dius at = 3.0 kpc for all primary galaxies; be- 
cause in the lMo et al.l (|l998l ) model, rd depends on the disc 
mass and on the halo spin parameter, we have choosen A 
in each case such that this value of is obtained in all 
primary systems. The corresponding spin parameters are 
-^0% = 0.033, A2o% = 0.034 and A4o% = 0.036. We consider 
two cases for the initial disc scale height: we adopt a "fidu- 
cial" value of zo = 0.4 kpc for direct comparison with P09, 
and also consider an initially thinner disc with zo = 0.25 kpc 
for two values of the gas fraction (0% and 20%). Edge-on 
surface brightness maps for the initial conditions of our pri- 
mary galaxies are shown in the upper and lower left panels 
of Figure [T] for zo — 0.4 and 0.25 kpc, respectively. 

The satellite systems consist only of a dark matter halo 
and a stellar bulge. The concentration parameter of the 
satellite halo is c = 11.98. Using the stellar-to-halo mass 



ratio we derive a stellar mass of Mb, sat ~ 6.3 x 10 Mq. We 
set the scale length of the bulge to rb,sat = 0.3 kpc. 

The primary systems always contain A'^dm = 4 x 10® 
dark matter, A^disc = 10® stellar disc and A^buigc = 5 x 10^ 
bulge particles. Models with 20% and 40% gas fractions 
contain A''gas,20% = 1-25 x 10^ and A''gas,40% = 3.33 x 10^ 
gas particles, respectively. The satellite system consists of 
A'dm = 9 X 10^ dark matter and A'^buigc = 10^ bulge parti- 
cles. We set the gravitational softening lengths to e = 50 pc, 
70 pc and 100 pc for stellar, gas and dark matter particles, 
respectively. 

Following P09 we choose orbits which are motivated 
by studies of substructure accretion in cosmological N- 
body simulations. The most likely values of the radial and 
tangential velocity components {vr and ft) are found to 
be respectively a t 90% and 60% of the virial velocity o f 
the primary halo (|Bensonll2005l : iKhochfar fc Burkertll2006l ') . 
For a halo of mass Mvir = 10^^ Mq the virial velocity is 
Wvir = 129kms~^ which results in an initial subhalo veloc- 
ity with Vr = 116kms~^ and vt = 77kms~^. The initial 
separation of the galaxies was chosen to be relatively large, 
dstart = 120 kpc, in order to prevent significant perturba- 
tions of the disc due to the sudden presence of the satel- 
lite's gravitational pull. Once dstart is chosen, the orbital 
parameters are uniquely determined. This, of course, limits 
our study to a single cosmologically motivated orbit. Our 
combination of initial velocity and distance fixes the peri- 
centric distance at ~ 18 kpc. We use a set of three orbital 
inclinations {9 = 60°, 45° and 30°, where 6 is the angle 
between the spin axes of the disc and the orbit) to investi- 
gate the effect of the inclination on the rate of thickening 
of the disc. All orbits are prograde and all simulations were 
evolved for a total of 6 Gyr. Prograde orbits are expected 
to be more destructive to the disc than retrograde mergers 
l| Velazquez fc WhitellT999f ). 

We summarize the parameters that are kept constant 
for all simulations in Table [1] and parameters that differ 
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Figure 1. Edge-on surface brightness maps for galaxies with an initial scale height of zq = 0.40 kpc (upper row) and zq — 0.25 kpc (lower row). The 
left column shows the initial models while the centre and the right columns show the final galaxies (t — 6 Gyr) for an initial gas fraction of 0% and 20%, 
respectively. A mass-to-light ratio of 3 (A-//L)q has been assumed, typical of the Milky Way in the B-band. 
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Figure 2. Density profiles for the initial (zq — 0.4 kpc) and final discs 
for an isolated galaxy and mergers with three inclination angles, for the 
simulations with no gas. The left panel shows the stellar mass density as 
a function of the distance to the galactic plane at a radius of 8 kpc. The 
right panel shows the cdgc-on projected surface brightness at a projected 
radius of 8 kpc. 

for the various simulation runs are summarized in Table [J] 
We label the diflterent simulations with the first letter I for 
isolated runs and an M for mergers. The second letter sig- 
nifies the gas fraction of the disc followed by a number sig- 
nifying the orbital inclination. We add a T for discs with 
an initial scale height of zq — 0.25 and an S for simula- 
tions with a lower star formation efficiency. We adopt a gas 
fraction-inclination combination of 20% and 60° as our fidu- 
cial model. 



3 RESULTS 

In order to study the evolution of the disc we compute the 
stellar mass density as a function of the distance to the 
galactic plane at a distance of Rq ~ 8 kpc from the disc 
spin axis. Figure [2] shows density profiles for a set of simu- 
lations that do not contain any gas. In the left panel, stellar 
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0.0 t , , . . , . . , , , . , . . , ] 

5 10 15 

R / kpc 

Figure 3. Disc scale height derived from fitting to the stellar mass density 
at different radii, for the initial (dotted line) and final states of a an isolated 
galaxy (lA; dashed line) and for a merger without gas (MA60, solid line) 
and with 20% gas {MB60; dot-dashed line), both for an inclination of 60°. 



mass density profiles of the initial disc (with zq = 0.4 kpc) 
and the final discs for an isolated galaxy and mergers with 
three inclination angles are plotted. To compare simulations 
to observations we also compute the edge-on projected sur- 
face brightness as a function of the distance to the galactic 
plane at a radius of Rq. This is done by taking a vertical 
slice at a distance of Rq to the galactic centre and assuming 
a mass-to-light ratio of M/L = 3(M/L)q, approp riate for 
the Milky Way in the B-band (|Zibetti et al.ll2009i ). In the 
right panel of Figure [2] the edge-on surface brightness pro- 
files are shown for the same simulations. As is clearly visible 
from both panels, the disc of the isolated galaxy is stable 
over the duration of the simulation, while during mergers 
the profiles clearly broaden, implying that the disc becomes 
thicker. 

We note that the projected edge-on surface density pro- 
files are thicker than the "slices" through the 3d mass den- 
sity profiles. To understand this we plot the disc scale heights 
derived by fitting to the 3d stellar mass density slices at dif- 
ferent radii in Figure [S] This shows that zo increases with 
increasing distance from the galaxy center. When computing 
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Figure 4. Evolution of the disc scale height for simulations with no gas. 
The solid line shows the scale height of an isolated galaxy while the other 
lines illustrate mergers with different inclinations. 



the projected edge-on surface density profiles, stellar parti- 
cles with a larger distance from the galactic centre than the 
specified radius are naturally included, since some of these 
particles happen to lie along the line of sight. These parti- 
cles are on average also more distant from the galactic plane, 
as Figure |3] indicates. This results in an apparently thicker 
edge-on surface density profile and a larger projected disc 
scale height. In order to compare our results with obser- 
vations we use the scale height obtained by fitting to the 
projected edge-on surface density profile in the following. 
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Figure 5. Evolution of the disc scale height for simulations with gas. 
The solid and the dashed lines show the isolated case while the other lines 
illustrate mergers with different gas fractions (0%, 20% and 40% for MA60, 
MB60 and MC60, respectively). 




3.1 Stability of the Initial Conditions and 
Evolution of the scale height 

Galactic discs are fragile systems and there is always the 
possibility that purely numerical effects can modify their 
morphology, e.g. through bar formation or flaring. The sta- 
bility of the initial disc is then a key point to be addressed 
before looking at the effects of satellite mergers. To quan- 
tify the amount of disc thickening we fit for the value of 
zo assuming a sech^ profile. We always fit to the projected 
edge-on surface brightness profile at Rq as discussed above. 

The dotted and dashed lines in Figure [3] show the val- 
ues of the disc scale height as a function of radius for an 
isolated galaxy (lA) ai t — and t — 6 Gyr, respectively. 
The final disc does not develop any appreciable flaring and 
the thickening at larger distances is negligible compared to 
the thickening due to satellite accretion events. In Figure [4] 
we show the resulting time evolution of zo for an isolated 
galaxy and mergers with varying inclinations, all contain- 
ing no gas. The isolated case (solid line) shows that in the 
absence of perturbations the thin disc is extremely stable 
during the 6 Gyr timescale over which we run the simula- 
tions. This is also true for the isolated runs including a gas 
component (IB and IC) as the solid and dashed lines in Fig- 
ure [S] show. This implies that any increase in zq is due to 
accretion events. We see that regardless of the inclination 
of the orbit, the disc thickens significantly in the merger 
simulations, with zq increasing by a factor of ~ 2. 

We note that we also ran a simulation of an isolated 
galaxy using one fourth the number of particles of our fidu- 
cial case (in all species) and find that the scale height in- 
creases by 50%. This demonstrates the importance of nu- 
merical resolution in determining the disc stability. 
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Figure fj. Evolution of the disc scale height for simulations without gas 
but with a more massive disc. The dotted lines shows the result for a sim- 
ulation with an additional collisionlcss disc having the same mass and dis- 
tribution as the gas disc in MC60. The dashed line shows the result of a 
simulation containing one stellar disc with increased mass (equal to the to- 
tal baryonic mass in MC60). The solid and dot-dashed lines indicate the 
reference cases (MA60, MC60, IC). We sec that the presence of dissipa- 
tional gas, not the more massive disc, is responsible for the reduction in 
disc thickening. 

3.2 The effect of gas in the disc 

All of the previous numerical experiments devoted to study- 
ing the heating and thickening of stellar discs by minor merg- 
ers have included only the dissipationless components of the 
galaxies: dark matter and stars. The presence of a gaseous 
component in the disc, which requires a hydrodynamical ap- 
proach to the problem instead of a purely gravitational one, 
has been neglected so far. The presence of gas may suppress 
disc thickening in two ways. One possibility is that the gas 
may absorb some of the kinetic impact energy of the merg- 
ing satellite, which can then be removed from the system 
by radiative cooling. In this way the impact energy that is 
transferred to the stars, causing heating and thickening, may 
be reduced. Another possible mechanism is that gas forms a 
new thin stellar disc after the merger. This disc could then 
cause the heated stars to contract towards the disc plane. 

We analyze the surface brightness profiles of the merger 
simulations with gas, and show the resulting evolution of 
the scale height in Figure [5] for initial gas fractions of 0%, 
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Figure 7. Star fo rmation rates for isolated galaxies and mergers, with 
initial gas fractions of 20% (IB, MB60, MB45, MB30) and 40% {IC, MC60) 
in the progenitor disc. 
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Figure 8. Evolution of the disc scale height for old stars (created in the 
initial conditions) and new stars (created during the simulation through 
star formation). The scale height for the combined sample is also shown for 
comparion. 



20% and 40%, and an inclination of 60° (MA60, MB60 and 
MC60). This shows that the presence of gas does indeed 
suppress the thickening of the disc by a minor merger. The 
final scale height increases by a factor of ~ 1.75 for the 20% 
gas case, and by only a factor of 1.5 for the 40% gas case, 
in contrast to the factor of 2 increase in the gas-free case 
(corresponding to a decrease in the final scale height of 25% 
and 50%, respectively). 

One might wonder whether the reduced thickening is 
due to the additional potential of the gas particles. In order 
to investigate whether the gas physics or the greater po- 
tential have a larger impact on the suppression of the disc 
thickening, we ran two additional coUisionless simulations. 
For the first simulation we use the initial conditions of MC60 
and convert all gas particles to stellar particles at t = Gyr 
(MC60X). We thus get an additional stellar disc which 
has the same profile and potential as the previous gaseous 
disc, but behaves coUisionlessly during the simulation. In 
a second simulation we create a galaxy like MA60 but in- 
crease the mass of the stellar disc to Mdisc = 4.0 x 1O^"M0 
(MA60Y). This galaxy has the same disc mass as MC60 
and MC60X. Figure [S] shows that in both coUisionless simu- 
lations (MC60X and MA60Y) the disc scale height increases 
by a factor of r^ 2. The final scale heights are equal or even 
larger than for the less massive galaxy without gas (MA60). 
This demonstrates that the suppression of the thickening in 
the runs including gas is not due to the additional potential 
of the gas disc but lies in the hydrodynamical nature of the 
particles. 

In order to gain insights into the physical process that 
is causing this change in behaviour, we examine the SFR as 
a function of time. The results are shown in Figure [7] The 
galaxy with a gas fraction of 40% (IC and MC60) naturally 
starts with a high SFR (~ 18 Mq yr~^), while the galaxy 
with a gas fraction of 20% (IB, MB60, MB45 and MB30) 
starts with a lower SFR of ~ 4 Mq yr~^. We see that 1:10 
mergers do not have a noticable impact on the SFR. 

The SFR quickly drops to low values under 1 Mq yr~^ 
(after 2.5 Gyr for a gas fraction of 40% and 1.5 Gyr for the 
gas fraction of 20%). This implies that most of the gas is con- 
sumed before the merger is complete. Thus, the galaxy is not 
able to reform a new thin stellar disc after the merger, which 
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Figure 9. Final edge -on projected surface brightness profile for the 40% 
gas case at a projected radius of 8 kpc for old, new and all stars. The solid 
line shows the fit to the combined stellar disc profile. 



could pull heated stars towards the galactic plane again. 
This means that the dominant process preventing disc heat- 
ing in the simulations above is the absorption of the kinetic 
impact energy of the satellite by the gas component. This 
can also be directly inferred from Figure \5\ from about 2 
Gyr, one can see that the growth of the scale height zq is 
slower in the simulations with gas, indicating that from the 
beginning the gas is absorbing the impact energy and pre- 
venting the disc from thickening considerably. If the main 
mechanism was the reforming of a new thin disc, one would 
see the disc thicken and then contract again. This would 
mean a different slope for the zo time evolution in runs with 
different gas fractions, however, we find similar slopes for all 
gas fractions indicating that this process is not dominant. 

Another way to demonstrate that the formation of a 
new disc does not reduce the disc thickening noticeably in 
our simulations is to compare the scale height of the stel- 
lar particles created in the initial conditions (old stars) to 
the scale height of the stellar particles that form during the 
simulation through SF (new stars). In Figure [S] we compare 
the evolution of zo for old and new stars in our fiducial case 
MB60. The new stars clearly form a thinner disc than the 
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Figure 10. Bvolut ion of the disc scale height for simulations with gas. 
The solid line shows the scale height of an isolated galaxy while the other 
lines illustrate mergers with different inclinations. 

old stars, however, the combined sample has a scale height 
which is only slightly thinner than the old stellar disc. This is 
due to the much larger mass contained in the old stellar disc 
component. Figure [5] shows the resulting edge-on projected 
surface brightness profile for the 40% gas case at 8 kpc for 
old, new and all stars and a fit to the combined stellar disc 
using a sech^ model. The combined profile does not show an 
obvious transition between the old and the new disc profile. 
However, one has to keep in mind that we are using a single 
sech^ model to fit two distinct profiles, implying that the 
fitted scale height of the combined sample decreases due to 
the presence of the new thin disc. The new and old discs 
have a final scale height of 0.38 and 0.66 kpc while the total 
stellar disc has zo — 0.61 kpc. This shows that even for a 
gas fraction of 40% (which is expected for disc galaxies at 
z ^ 1 but not for local galaxies at the present time), the new 
thin disc is not massive enough to reduce the total disc scale 
height significantly. In order to form a new thin disc with a 
mass comparable to the old stellar disc, the SFR after the 
merger would have to be much higher. This could occur in 
the case of an extremely g as-rich initial d isc 90% gas; 
[Robertson eT al. 2006; Hop kins et al]|2008l '). 

In order to investigate the influence of the orbital incli- 
nation on the disc thickening we compute the evolution of 
the scale height for mergers with gas fractions of 20% for 
three different orbital inclinations 6 = 60° , 45°, 30° (MB60, 
MB45, MB30 respectively). We plot the results in Figure [TOl 
The black solid line (IB) represents the isolated case with 
gas (20%), which is seen to be very stable as in the no-gas 
case. The results for different orbital inclinations show that 
the value of has a small effect on the growth of the scale 
height. Regardless of the inclination, the overall increase of 
the scale height for simulations including gas is about 25% 
smaller than in the respective dissipationless simulations. 

3.3 A thinner initial disc 

The choice of initial scale height for our primary disc was 
somewhat arbitrary. In order to investigate whether a thin- 
ner initial disc would result in a thinner disc today, we 
simulate a merger of a disc with an initial scale height of 
zo = 0.25 kpc both without and with gas (20%) and an in- 
clination of 60°. The resulting evolution of zq is shown in 
Figure 1111 The final scale height of the disc without gas is 
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Figure IX. Evolution of the scale height for a disc with a smaller initial 
scale height. The solid line shows the 'thin disc' simulation without gas, the 
dotted line shows thin disc case with 20% gas, and the dashed line shows 
the 20% gas case with the original choice of initial scale height. 

as large as the scale height of the corresponding initially 
thicker discs. This may be due to the fact that an initially 
thicker disc is more rob ust to heating by accretion events 
(IKazantzidis et al .12009I ). When gas is present, the final scale 
height is slightly smaller in the case with an initially thin- 
ner disc, but even a thinner initial disc is transformed into 
a system with a scale height of more than 0.6 kpc. Thus, 
even when the effects of gas in the progenitor disc are in- 
cluded, it seems to be difficult to obtain a disc with a scale 
height that is substantially smaller than this value following 
a minor (mass ratio greater than 1:10) merger. 

3.4 Lower SF efficiency 

In order to investigate the sensitivity of our results to the pa- 
rameters chosen for star formation and SN feedback, we have 
run some simulations with lower star formation efficiency. 
In particular, we were interested to see whether adopting a 
lower star formation efficiency would allow a more significant 
new thin disc to reform following the merger. For this test 
we employ a star formation timescale = 8.4 Gyr, a cloud 
evaporation parameter Ao — 4000 and a SN "temperature" 
TsN = 4.0 X 10* K in order to h ave a SFR of ^ 1 Ma yr'^ 
for a Milky Way-hke Galaxy l|Springel et al.1 bOOsh . We 
run this case for initial scale heights of zo = 0.4 kpc and 
Zo = 0.25 kpc, in both cases using our fiducial value of a 
20% gas fraction and an inclination of 60°. 

We first compare the SFRs of the simulations with low- 
ered SF efficiencies with our fiducial results. The left panel 
of Figure [12] shows the SFRs for the simulations MB60, 
MB60T, MB60S and MB60TS. As we see, the SFRs in 
MB60S and MB60TS are initially lower than those of MB60 
and MB60T, respectively. After the merger {t ~ 1.8 Gyr) 
however, the simulations with the lower SF efficiency have a 
higher SFR, since MB60 and MB60T have already consumed 
most of their gas supply. Due to the lower SFR during the 
early stages of the simulation, the amount of newly formed 
stellar mass is still lower in MB60S and MB60TS, as can 
be seen in the right panel of Figure 1121 This means that 
the thin disc consisting of new stars affects the overall disc 
scale height even less than in the simulations MB60 and 
MB60T. The resulting scale heights are plotted in Figure 
[T3l As expected from the SFRs, the scale heights of MB60S 
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Figure 12. Left panel: SFR tor simulations with the default (MB60, 
MB60T) and lowered SF efficiencies (MB60S, MB60TS), for the disc with 
initial scale height 20 — ^^'4 kpc and zq — 0.25 kpc (T). Right panel: Mass 
of stellar particles formed during the simulations. 
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Figure 13. Evolut ion of the disc scale height for simulations with higher 
{MB60, MB60T) and lower SF efficiencies {MB60S, MB60TS). 

and MB60TS are slightly larger than those of MB60 and 
MB60T, respectively. 

Clearly, in order for the disc thickening problem to be 
solved by the growth of a new stellar disc following the 
merger, the mass of the new disc must be comparable in 
mass to that of the old stellar disc. This implies that a sig- 
nificant amount of gas must remain available at the end of 
the merger. Based on our results, when the SF and SN feed- 
back parameters are tuned be consistent with the empirical 
Kennicutt law, most of the gas tends to be consumed during 
the merger. Simulations with much higher initial gas frac- 
tions (~ 90%) or which included a new source of cold gas 
(e.g. via cooling and accretion from a hot halo) might be 
able to regrow a more massive thin disc. 



4 CONCLUSIONS AND DISCUSSION 

We investigated the role of dissipational gas physics in the 
vertical heating and thickening of disc galaxies by minor 
mergers (mass ratio 1:10). We used the parallel TreeSPH- 
code GADGET-2 to simulate a suite of minor merger simu- 
lations for a Milky Way-like primary galaxy. The suite con- 
sists of coUisionless simulations as well as runs containing 
gas in the disc. Using cosmologically motivated orbital pa- 
rameters, we ran simulations with three different orbital in- 



clinations (60°, 45° and 30°) and two different initial disc 
scale heights {zo = 0.4 kpc and zo — 0.25 kpc). 

We fit for the value of the scale height zo, assuming a 
sech^ form for the vertical density profile. We showed that 
the scale height derived from a fit to the mass density is a 
strong function of the radial distance to the disc axis: as is 
well known, minor mergers tend to cause flaring in the outer 
part of the disc. This results in projected edge-on surface 
density profiles which appear thicker than the actual three 
dimensional density profiles since stars that are at a large 
distance from the galactic centre and thus also more distant 
from the galactic plane are viewed along the line of sight. 

We found that in dissipationless simulations, minor 
mergers (without gas) caused the scale height of the disc 
to increase by about a factor of two, with the orbital incli- 
nation affecting the amount of thickening by only about five 
percent. Thus, in qualitative agreement with the results of 
POO, we find that in the absence of gas, cold thin discs are 
destroyed by minor mergers. However, it is interesting that 
in spite of the fact that the parameters of our simulations 
were chosen to be very similar to those of P09, they found 
that the scale height of the disc increased by a larger factor 
of ~ 3, leading to a final disc about twice as thick as ours 
(despite having the same initial scale height). 

Although the discs in both simulation suites are stable 
against secular evolution effects, this does not necessarily 
ensure that both models are equally resistant to strong per- 
turbative effects during an accretion event. Possible reasons 
for this discrepancy may be the use of a different IC genera- 
tor and a different simulation code. The velocity dispersions 
for example are fixed differently; while in our ICs the radial 
and the vertical velocity dispersions are equal, the code used 
by POO reproduces the velocity ellipsoid for the Milky Way. 
This may result in an initially hotter disc for our typical 
galaxy (for a Mvir = 10^^ Mq halo) which is more resilient 
to massive accretion events. Other intrinsic differences in 
the simulations of POO are a more massive disc (which may 
potentially be less stable due the formation of a bar) and a 
larger stellar mass of the satellite (which results in a larger 
gravitational pull each time the subhalo passes the disc) . We 
also note that POO use a different concentration parameter 
(for both the primary halo and the subhalo) and chose a 
Sersic profile for the bulge and an NFW profile for the halo 
component while we use a Hernquist profile for both compo- 
nents. It is not obvious which of the mentioned differences 
has the largest impact on the resulting disc scale heights. 

We note that our re sults are consistent with those of 
IViUalobos fc Helmil | 20081 ). These authors start with similar 
initial conditions but employ a thicker disc {zq = 0.7 kpc). 
The final scale heights are between 1.0 and 1.3 kpc, depend- 
ing on the orbital inclination and are not as extreme as found 
by POO. The relative thickening is even lower than in our 
simulations, confirming that initially thicker discs are more 
robust to accretion events. 

We also investigated mergers in which the progenitor 
disc initially contained a gas fraction of 20 or 40 %. We 
found that the presence of gas reduces the final scale heights 
by 25% (50%) for a gas fraction of 20% (40%). The final 
scale heights were between 0.6 and 0.7 kpc (1.5 to 1.75 times 
larger than the initial value), depending on the initial gas 
fraction and the orbital inclination. We argued that the pres- 
ence of gas can have an impact on disc thickening via two 
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different mechanisms. One process is the absorption of ki- 
netic impact energy by the gas. This energy can then be 
dissipated via radiation. Another possible effect is the for- 
mation of a new thin disc, which can cause heated stars to 
recontract towards the disc plane. We showed that in our 
simulations, in which the SF and SN feedback parameters 
were set to be consistent with the empirical Kennicutt law 
for the initial disc, most of the gas is consumed during the 
merger, and therefore the regrowth of new thin discs has a 
negligible impact on the mass-weighted scale height of the 
post-merger galaxy. Therefore, it seems that the main pro- 
cess that suppresses disc thickening in the presence of gas 
in our simulations is the absorption of impact energy by the 
gas. We intend to investigate the details of the physics of 
this process in more detail in future work. 

We computed scale heights for old stars (i.e. stellar par- 
ticles that were present in the initial conditions) and for new 
stars (stellar particles created during the simulation through 
SF). The scale height for old stars was found to be slightly 
higher than the overall value (~ 0.73 kpc vs ~ 0.67 for our 
fiducial case of 20% gas and a 60° inclination) while the 
scale height of the new stars was about half of that value 
(~ 0.35 kpc). We showed that the final mass of the new 
disc is small compared to that of the old disc. This indicates 
that although a new thin disc does form, it does not have a 
significant effect on the final thickness of the total disc. We 
argue that in order to have a noticeable effect, the new thin 
disc would need to be comparable in mass to the old disc, 
which would require that a significant mass of cold gas is 
still present at the end of the merger. 

We ran simulations with two different initial scale 
heights, zo = 0.40 kpc (fiducial case) and zo = 0.25 kpc 
(thin). As it turns out, thinner discs are more unstable to 
heating and are therefore thickened more by the merger than 
initially thicker discs. As a result, the two cases result in 
discs with nearly the same final scale height (~ 0.6 kpc for 
the 20% gas fiducial case). 

To study the sensitivity of our results to the param- 
eters controlling SF and SN feedback, we ran simulations 
with a lower star formation efficiency. In these simulations, 
less gas is consumed and so the SFR at the end of the 
merger is higher. However, we found that the final mass 
of the new disc was still lower than in the simulations with 
the higher SF efficiency. This results in final scale heights 
that are slightly larger than in the fiducial case. We con- 
clude that in order to reform a new thin disc comparable in 
mass to the old disc, either the initial gas fraction would 
have to be much higher, or an external fueling reservoir 
(such as cooling and accretion from a hot halo) would be 
needed (ISommer-Larsen et al. 20031: Kaufmann et al.|[2006l : 
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In light of our results, we can now reassess whether 
disc thickening by minor mergers presents a serious prob- 
lem for CDM. First, we need to compare the fraction of 
disc galaxies that are expected to have had a minor merger 
(mass ratio greater than 1:10 but less than about 1:4) since 
2: ~ 1 with the fraction of observed galaxies with thin 
discs. Observations show that roug hly 70% of Milky Way- 



1:10 or larger in the past 8 Gyr. This implies that 30% of all 
Milky Way-like galaxies did not have such a merger, and as a 
result can be expected to have retained their thin disc. This 
means that if less than 30% of all Milky Way-like galaxies 
are found to have scale heights that are substantially lower 
than those found in our simulations (~ 0.6 kpc) then there 
is no discrepancy. 

Obtaining statistically unbiased observational measure- 
ments of scale heights for a complete sample of galax- 
ies is difficult, due to small sample sizes, dust extinc- 
tion, and inclination effects. Still there have been many 
studies focusing on the vertical structure of galaxies (e.g . 
Shaw fc Gilmor c 1989. 1990:' de Griis fc van der Kruitll99e : 
de Griis et anil997: .Pohlen et all l2000l: iKregel et all I2OO2 



Bizvaev fc KaisinI 20041 ). All these studies suggest that the 



majority of observed galaxies with properties similar to 
those of the Milky Way have sech^ scale heights in the range 
of 0.6 < 20 < 1.0 kpc with a ratio between scale height and 
scale length between 0.2 and 0.3. These values have been 
determined from observed two-dimensional surface bright- 
ness profiles. Usually it is assumed that scale length and 
scale height are independent parameters. A model 2d sur- 
face brightness profile consisting of an exponential radial 
factor and a vertical factor (sech^, sech or exp) is then fit- 
ted to the observed 2d profile. All scale heights cited here 
have been converted to a sech^ profile. 

A statistical study o f the vertical structure of sp i- 
ral galaxies is presented in [Schwarzkopf fc Dettrriail l|2000l ). 
They found that the majority (~ 60%) of galaxies have verti- 
cal scale heights less than 1.1 kpc, with a maximum between 
0.4 kpc zo =^ 0.8 kpc, while galaxies with a very thin disc 
(20 < 0.4 kpc) are extremely rar e . Sim ilar results were ob- 
tained bv lYoachim fc DalcantonI (|2006l ). who derived scale 
heights for a sample of edge-on galaxies and presented a rela- 
tion between scale height and circular velocity. They showed 
that Galaxy-sized systems are expected to have scale heights 
in the range 0.6 kpc ^ 20 ^ 1.2 kpc. In order to compare our 
results to these studies, we compute 2d edge-on disc surface 
brightness profiles (as shown in Figure [1]) and fit a model 2d 
profile consisting of an exponential radial and sech^ vertical 
disc to the simulated profiles. For the coUisionless fiducial 
run MA60 the final scale height is 0.74 kpc, while for the 
fiducial runs including a gas disc MB60 and MC60 the scale 
heights are 0.66 and 0.62 kpc, respectively. These values fall 
precisely in the observed range. 

Previous theoretical studies, such as POO, have fo- 
cussed on whether the observed scale height of the Milky 
Way is in confiict with the predictions of CDM models. 
We revisit this issue as well. There have been various 
studies of the vertical structure of the Milky Way (e.g. 
Bahcall fc Soneiralll980l:lKent et al]ll99ll:lRerd fc Maiewskil 
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studies find exponential scale heights for the Milky Way's 
old thin disc of ~ 0.3 kpc and for the young star forming 
disc of /iz ~ 0.1 kpc. As is standard in the literature, we con- 
vert these exponential scale heights to values corresponding 
to a sech^ function using 20 = 2/iz (which yields very good 
agreement between the two profiles at 2 J5 20 and reason- 
able agreement at zo > 2). This results in scale heights of 
zo = 0.6 kpc (old disc) and zo = 0.2 kpc (young disc). The 
former is just what we find for the overall scale height in 
our fiducial merger with gas, while the latter is very close to 
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what we find for the disc of new stars, as shown in Figure [S] 
Thus, in hght of our results, even if the Milky Way did expe- 
rience a recent 1:10 merger, if the disc contained even 20% 
gas (similar to the gas fraction today) at the time of the 
merger, there is no conflict with observations. Of course, it 
is also possible to explain the Milky Way by assuming that 
it is one of the 30% of systems that did not e xperience a 1:10 
merge r. It has been previously suggested bv lHammer et ahl 
((20071), based on the unusually low mass and angular mo- 
mentum of the Milky Way disc, that our Galaxy may have 
had a particularly quiescent formation history. 

We also note that while we are fitting the scale height to 
mass profiles, observational studies use luminosity profiles 
to der ive zo . Based on the models of iBruzual Sz CharlotI 
l|2003t ) for a single age stellar population, we expect the B- 
band mass-to-light ratio of very young stars {ta < 2 Gyr) 
and old stars {ta > 8 Gyr) to differ by a factor of ~ 4. 
This implies that the new stellar particles should carry a 
correspondingly larger weight in a luminosity- weighted fit, 
resulting in a smaller estimated scale height. However, as 
seen in Figure 1121 the mass of such young stars formed in 
our simulations is very small (only about a tenth of the mass 
of the pre-existing stellar disc), so even if we accounted for 
the age-dependent mass-to-light ratio, the overall measured 
scale height would not differ significantly from the mass- 
weighted one. Moreover, young stars will also tend to be 
more enshrouded by dust, further reducing their impact on 
the measured scale height. 

Finally, the simulations presented here still neglect the 
larger scale environment and cosmological growth of the 
galaxy. The progenitors of Milky Way discs presumably had 
somewhat different properties (smaller masses, higher gas 
fractions, and possibly thicker discs) than their present-day 
counterparts. Perhaps most importantly, one expects these 
systems to have accreted a significant amount of new mate- 
rial over the course of the past 6-8 Gyr via accretion and 
cooling from a hot gas halo. This additional supply of cold 
gas to the disc could enhance the formation of a new, thin 
disc and further reduce the thickening by minor mergers. We 
plan to re-examine this problem in the context of simula- 
tions with a cosmologically motivated merger and accretion 
history, and including cooling from a hot halo, in a future 
work. 

In final summary, in contrast with the conclusions 
reached by P09 using dissipationless simulations, we con- 
clude that the existence of a thin disc in the Milky Way 
is not in obvious confiict with the predictions of the CDM 
model. Although the Milky Way may have had an unusually 
quiescent merger history, mooting the issue of disc thicken- 
ing by minor mergers, our results indicate that one does 
not need to rely on such an argument. When we include a 
moderate amount of cold gas in the progenitor disc (similar 
to the gas fractions observed in typical spirals today), the 
scale heights of simulated discs that have experienced a 1:10 
merger are in good agreement with the observed values both 
for the Milky Way thin disc and for external galaxies. 
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